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Abstract 



The mean value locates the center of the probability mass distribution induced by X on the real line. 
In this unit, we examine how expectation may be used for further characterization of the distribution 
for X. In particular, we deal with the variance and its square root the standard deviation. We identify 
some important properties of variance and identify the concept covariance. The variance is calculated 
for several distributions, comparing analytical and Matlab results. 

In the trcatmont of tho mathematical expection of a real random variable X, we note that the mean 
value locates the center of the probability mass distribution induced by X on the real line. In this unit, we 
examine how expectation may be used for further characterization of the distribution for X. In particular, we 
deal with the concept of variance and its square root the standard deviation. In subsequent units, we show 
how it may be used to characterize the distribution for a pair {X, Y} considered jointly with the concepts 
covariance, and linear regression 

1 Variance 

Location of the center of mass for a distribution is important, but provides limited information. Two markedly 
different random variables may have the same mean value. It would be helpful to have a measure of the 
spread of the probability mass about the mean. Among the possibilities, the variance and its square root, 
the standard deviation, have been found particularly useful. 

Definition. The variance of a random variable X is the mean square of its variation about the mean 
value: 



The standard deviation for X is the positive square root ax of the variance. 
Remarks 

• If X {(jj) is the observed value of X, its variation from the mean is X (w) — ^x- The variance is the 
probability weighted average of the square of these variations. 

• The square of the error treats positive and negative variations alike, and it weights large variations 
more heavily than smaller ones. 

• As in the case of mean value, the variance is a property of the distribution, rather than of the random 
variable. 

• We show below that the standard deviation is a "natural" measure of the variation from the mean. 
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Var [X] = a\=E {X - fixf where nx = E [X] 
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• In the treatment of mathematical expectation, we show that 



E 



{X - cf is a minimum iff c=E[X], in which case E {X - E [X]f = E [X"^] -E"^ [X] (2) 



This shows that the mean value is the constant which best approximates the random variable, in the 
mean square sense. 

Basic patterns for variance 

Since variance is the expectation of a function of the random variable X, we utilize properties of expec- 
tation in computations. In addition, we find it expedient to identify several patterns for variance which are 
frequently useful in performing calculations. For one thing, while the variance is defined as {X — jix)^ , 
this is usually not the most convenient form for computation. The result quoted above gives an alternate 
expression. 

(VI): Calculating formula. Var [X] = E [X^] - E^ [X]. 

(V2): Shift property. Var [X + 6] = Var [X]. Adding a constant b to X shifts the distribution (hence its 
center of mass) by that amount. The variation of the shifted distribution about the shifted center of 
mass is the same as the variation of the original, unshifted distribution about the original center of 
mass. 

(V3): Change of scale. Var [aX] = a^Var [X]. Multiplication of X by constant a changes the scale by a 
factor \a\. The squares of the variations are multiplied by a^. So also is the mean of the squares of the 
variations. 

(V4): Linear combinations 

a. Var [aX ± bY] = a^Var [X] + b^Yar [Y] ± 2ab {E [XY] -E[X]E [Y]) 

b. More generally. 



Var 



"^akXk 



.fc=i 



J2 4 Var [Xk] + 2 ^ a,aj {E [X^Xj] - E [Xi] E [Xj]) 



(3) 



fc=i 



The term Cjj = E [XiXj] — E [Xi] E [Xj] is the covariance of the pair {X^, Xj}, whose role we study 
in the unit on that topic. If the Cij are all zero, we say the class is uncorrelated. 

Remarks 

• If the pair {X, Y} is independent, it is uncorrelated. The converse is not true, as examples in the next 
section show. 

• If the ttj = ±1 and all pairs are uncorrelated, then 



Var 



y^^ajXi 



.k=l 



^Vax[X, 



(4) 



fe=i 



The variance add even if the coefficients are negative. 



We calculate variances for some common distributions. Some details are omitted — usually details of algebraic 
manipulation or the straightforward evaluation of integrals. In some cases we use well known sums of infinite 
series or values of definite integrals. A number of pertinent facts are summarized in Appendix B^. Some 
Mathematical Aids. The results below are included in the table in Appendix C^. 
Variances of some discrete distributions 

^"Appendix B to Applied Probability: some mathematical aids" <http://cnx.oig/content/m23990/latest/> 
^"Appendix C: Data on some common distributions" <http://cnx.org/content/m23992/latest/> 
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1. Indicator function X = IeP {E) =p,q=l— pE [X] = p 

E[X^]-E^X] = E[ll] -p^ = E[lE]-p^=p-p^=p{l-p)=pq (5) 

2. Simple random variableX = Y17=i^i-^^i (primitive form) P (Ai) = pi- 

n 

Var [X] = ^ tiPiQi - 2 ^ titjPiPj,sinceE [/aJaJ = Oi ^ j (6) 

i=l i<j 

3. Binoinial{n,p). X = 7B,with{/B, : 1 < i < njiidP {Ei) = p 

n n 

Var [X] = [^^J = ^ = (7) 



i=l 



4. Geometnc{p). P {X = k) = pq'^^k > OE [X] = q/p 
We use a trick: [X^] = ^ [X (X - 1)] + [X] 



OO CO O 2 

; n I o (-^ ^ 



/s=0 fc=2 



Var [X] = 2^ + q/p - {q/pf = q/p^ (9) 

5. Poisson{iJ,)P {X = k)= e"'^ ^Vfc > 

Using E [X2] = E[X {X - 1)] + E[X], Mve have 

oc ^ CO k—2 

fc=2 ■ k=2 ^ '' 

Thus, Var [X] = /i^ + /x — /z^ = /x. Note that both the mean and the variance have common value ji. 
Some absolutely continuous distributions 
1. Uniform on (a, b)fx{t) = ^a<t<bE [X] = 



2 



E m = / t^dt = ° so Var [X] 

L J 3(6-a) ^ ^ 



3 (6 - a) 4 12 



(11) 



2. Symmetric triangular {a, b) Because of the shift property (V2) ("(V2)", p. 2), we may center the 
distribution at the origin. Then the distribution is symmetric triangular (— c, c), where c = (6 — a) /2. 
Because of the symmetry 

Var [X] = E [X^] = J' t'fx {t) dt = 'if^ ffx {t) dt (12) 



Now, in this case. 



{t) = ^0 < t < cso thatE [X2] = ^ {ct^ -t^)dt=j= (13) 
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3. Exponential (A) fx {t) = Xe-^*,t> QE[X] = 1/A 

E [X^] = / Xfe-^^dt = — sothatVax [X] = l/A^ (14) 
Jo X 

4. Gamma{a, X)fx (<) = j^A^t^-^e"^** > OE [X] = f 

Hence Var [X] = a/X^. 

5. NoTmal{n,a'^)E[X] = 

Consider Y ^ N ,E [Y] = 0, Var [Y] = ^ t^e'^^/'^dt = 1. 

X = aY + /liimpliesVar [X] = a^Var [Y] = cr^ (16) 

Extensions of some previous examples 

In the unit on expectations, we calculate the mean for a variety of cases. We revisit some of those 
examples and calculate the variances. 

Example 1: Expected winnings (Example 8^ from "Mathematical Expectation: Simple 
Random Variables") 

A bettor places three bets at $2.00 each. The first pays $10.00 with probability 0.15, the second 
$8.00 with probability 0.20, and the third $20.00 with probability 0.10. 
SOLUTION 

The net gain may be expressed 

X = 10/4 + 8/b + 20/c - 6, with P (A) = 0.15. P (B) = 0.20, P (C) = 0.10 (17) 

We may reasonbly suppose the class {A, B, C} is independent (this assumption is not necessary 
in computing the mean) . Then 

Var [X] = lO^P (A) [1 - P (A)] + S^P (B) [1 - P (B)] + 20^ P (C) [1 - P (C)] (18) 
Calculation is straightforward. We may use MATLAB to perform the arithmetic. 

c = [10 8 20] ; 
p = 0.01*[15 20 10] ; 

q = 1 - p; 

VX = sum(c.~2.+p.+q) 
VX = 58.9900 

Example 2: A function of X (Example 9^ from "Mathematical Expectation: Simple 
Random Variables") 

Suppose X in a primitive form is 

X = -3Ic, - + '^Ics - 3Ic, + "^Ic, - Ic, + Ic, + 'ilcs + 3Ic, + 2/cio (19) 

with probabilities P{Ci) = 0.08, 0.11, 0.06, 0.13, 0.05, 0.08, 0.12, 0.07, 0.14, 0.16. 

Let g {t) =t^ + 2t. Determine E [g {X)] and Var [g {X)] 



^"Mathematical Expectation: Simple Random Variables", Example 8: Expected winnings 
<http://cnx.Org/content/m23387/latest/#fs-idll64693988060> 

^"Mathematical Expectation: Simple Random Variables", Example 9: Expectation of a function of X 
<:http://cnx.org/content/m23387/latest/#fs-id2920427> 
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c = [-3 -12-3 4 -1 


1 2 3 2]; 


7. Original coefficients 


pc = 0.01*[8 11 6 13 5 8 


12 7 14 16] ; 


y, Probabilities for C_j 


G = c.~2 + 2*c 




'/. e(c i) 


EG = G*pc' 




°/o Direct calculation E[g(X)] 


EG = 6.4200 






VG = (G.-2)*pc' - EG-2 




°/o Direct calculation Var[g(X)] 


VG = 40.8036 






[Z,PZ] = csort(G,pc) ; 




y, Distribution for Z = g(X) 


EZ = Z*PZ' 




y. E[Z] 


EZ = 6.4200 






VZ = (Z.~2)*PZ' - EZ~2 




y. Var[Z] 


VZ = 40.8036 







Example 3: Z = g (X, Y) (Example 10^ from "Mathematical Expectation: Simple Ran- 
dom Variables") 

We use the same joint distribution as for Example 10^ from "Mathematical Expectation: Simple 
Random Variables" and let g (t, u) = t'^ + 2tu — 3u. To set up for calculations, we use jcalc. 

jdemol y. Call for data 

jcalc y. Set up 

Enter JOINT PROBABILITIES (as on the plane) P 
Enter row matrix of VALUES of X X 
Enter row matrix of VALUES of Y Y 

Use array operations on matrices X, Y, PX, PY, t, u, and P 
G = t."2 + 2+t.+u - 3+u; °A Calculation of matrix of [g(t_i, u_j)] 
EG = total(G.*P) y. Direct calculation of E[g(X,Y)] 

EG = 3.2529 

VG = total(G.-2.+P) - EG-2 I Direct calculation of Var[g(X,Y)] 
VG = 80.2133 

[Z,PZ] = csort(G,P); °A Determination of distribution for Z 

EZ = Z*PZ' y, E[Z] from distribution 

EZ = 3.2529 

VZ = (Z.-2)+PZ' - EZ~2 y. Var[Z] from distribution 

VZ = 80.2133 

Example 4: A function with compound definition (Example 12^ from "Mathematical 
Expectation; General Random Variables") 

Suppose X ^ exponential (0.3). Let 

for X < 4 

Z = { - =7[o,4](X)X' + /(4,cx>](^)16 (20) 

16 for X > 4 ^ ' ^ ' 

Determine E[Z] and Var [Z]. 
ANALYTIC SOLUTION 

E [g (X)] = Jg{t) fx {t) dt = /[o,4] (t) f0.3e-°''' dt + 16^; [/(4,oo] {X)] (21) 



^"Mathematical Expectation: Simple Random Variables", Example 10: Expectation for Z = p {X, Y) 
<http://cnx.Org/content/m23387/latest/#fs-idll64693552188> 

^"Mathematical Expectation: Simple Random Variables", Example 10: Expectation for Z = p {X, Y) 
<http://cnx.Org/content/m23387/latest/#fs-idll64693552188> 

'^"Mathematical Expectation; General Random Variables", Example 12: A function with a compound definition: truncated 
exponential ^http://cnx.org/content/m23412/latest/#fs-id6697312> 
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= / i^O.Se""-^* dt + 16P (X > 4) w 7.4972 (by Maple) (22) 
Jo 

= J[o,4] {X) + /(4,oo] {X) 256 (23) 

256E [1(^4,00] {X)] = / t^0.3e-°-^* dt + 256e-^ '^ « 100.0562 (24) 
Jo 



/[o,4] {t) f''0.3e-°-2* 

Var [Z] = E [Z'^] - E"^ [Z] w 43.8486 (by Maple) 



(25) 



APPROXIMATION 

To obtain a simple aproximation, we must approximate by a bounded random variable. Since 
P (X > 50) = e~i^ « 3 • 10"'^ we may safely truncate X at 50. 

tappr 

Enter matrix [a b] of x-range endpoints [0 50] 
Enter number of x approximation points 1000 
Enter density as a function of t . 3+exp(-0 . 3+t) 
Use row matrices X and PX as in the simple case 

M = X <= 4; 

G = M.+X.-2 + 16+(1 - M): 
EG = G*PX' 
EG = 7.4972 
VG = (G.-2)+PX' - EG-2 
VG = 43.8472 
[Z,PZ] = csort(G,PX) ; 
EZ = Z*PZ' 
EZ = 7.4972 
VZ = (Z.-2)+PZ' - EZ~2 
VZ = 43.8472 



g(X) 
E[g(X)] 

7. Var[g(X)] 

'/. Theoretical = 43.8486 

Distribution for Z = g(X) 
°/o E[Z] from distribution 

Var[Z] 



Example 5: Stocking for random demand (Example 13^ from "Mathematical Expec- 
tation; General Random Variables") 

The manager of a department store is planning for the holiday season. A certain item costs c 
dollars per unit and sells for p dollars por unit. If tlio domand exceeds the amount m ordered, 
additional units can be special ordered for s dollars per unit (s > c). If demand is less than the 
amount ordered, the remaining stock can be returned (or otherwise disposed of) at r dollars per 
unit (r < c). Demand D for the season is assumed to be a random variable with Poisson (/i) 
distribution. Suppose /x = 50, c = 30, p = 50, s = 40, r = 20. What amount m should the manager 
order to maximize the expected profit? 
PROBLEM FORMULATION 

Suppose D is the demand and X is the profit. Then 

For D < m, X = D (p — c) — (m — D) [c — r) — D {p — r) + m{r — c) 
For D > m, X = m{p — c) + {D ~ m) {p ~ s) = D {p — s) + m {s — c) 

It is convenient to write the expression for X in terms of Im, where M = (— oo,m]. Thus 

X = Im {D) [D{p-r) + m{r- c)] + [1-Im {D)] [D {p - s) + m{s - c)] (26) 

® "Mathematical Expectation; General Random Variables", Example 13: Stocking for random demand (see Exercise 4^ 
from "Problems on Functions of Random Variables") <http://cnx.Org/content/m23412/latest/#fs-idl0517737> 
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Then 



= D {p — s) + m{s — c) + Im {D) [D [p — r) + m{r — c) — D [p — s) — m{s — c)] 

= D {p - s) + m{s - c) + Im {D) (s -r)[D- m] 

E [X] = ip - s) E[D] + m{s - c) + {s - r) E [Im {D) D] - {s - r) mE [Im {D)] 



(27) 
(28) 

(29) 



We use the discrete approximation. 
APPROXIMATION 



> mu = 50; 

> n = 100; 

> t = 0:n; 

^ pD = ipoisson(mu,t) ; 
30; 



y. Approximate distribution for D 



> c 

> P 

> s 

> r 
» m 
» for i 



55; 

1 : length (m) 



M = t<=m(i) ; 

G(i,:) = (p-s)*t + m(i)*(s-c) + (s-r)*M.*(t 



7, Step by step calculation for various m 
m(i)); 



end 



> 


EG 


= G*pD 


5 . 

J 










> 


VG 


= (G.- 


2)+pD' - 


EG 


-2; 




SG 


=sqrt (VG) ; 










> disp([EG' 


;VG' 


;SG'] 


') 








l.Oe+04 * 

















0931 


1. 


1561 







0108 







0936 


1. 


3117 







0115 







0939 


1. 


4869 







0122 







0942 


1. 


6799 







0130 







0943 


1. 


8880 







0137 







0944 


2. 


1075 







0145 







0943 


2. 


3343 







0153 







0941 


2. 


5637 







0160 







0938 


2. 


7908 







0167 







0934 


3. 


0112 







0174 







0929 


3. 


2206 







0179 



Example 6: A jointly distributed pair (Example 14^'' from "Mathematical Expectation; 
General Random Variables") 

Suppose the pair {X, Y} has joint density fxY {t, u) = 3u on the triangular region bounded by 
u = 0, u=l + t, u=l-t. Let Z = g{X,Y) = X^ + 2XY. 
Determine E [Z] and Var [Z] . 
ANALYTIC SOLUTION 



^''''Mathematical Expectation; General Random Variables", Example 14: A jointly distributed pair 
<http://cnx.Org/content/m23412/latest/#fs-id7912879> 
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E[Z] = / / {t^ + 2tu) fxY (t, u) dudt = 3 u {e + 2tu) dudt + 

3 /o "* u {t^ + 2iii) dudt = 1 /lO 



/O rl+t pi pl-t 

/ + 2tu)^dudt + 3 / / u{t'^ + 2tu)'^ dudt = 3/35 

-1 Jo J a J a 

Vax [Z] = E [Z^] - E^ [Z] = 53/700 w 0.0757 



APPROXIMATION 
tuappr 

Enter matrix [a b] of X-range endpoints [-1 1] 

Enter matrix [c d] of Y-range endpoints [0 1] 

Enter number of X approximation points 400 

Enter number of Y approximation points 200 

Enter expression for joint density 3*u. + (u<=min(l+t , 1-t) ) 

Use array operations on X, Y, PX, PY, t, u, and P 

G = t.-2 + 2*t.+u; 7. g(X,Y) = X-2 + 2XY 

EG = total (G.*P) E[g(X,Y)] 

EG = 0.1006 '/. Theoretical value = 1/10 

VG = total(G.-2.*P) - EG~2 

VG = 0.0765 Theoretical value 53/700 = 0.0757 

[Z,PZ] = csort(G,P); Distribution for Z 

EZ = Z*PZ' E[Z] from distribution 

EZ = 0.1006 

VZ = Z.~2*PZ' - EZ~2 

VZ = 0.0765 

Example 7: A function with compound definition (Example 15^^ from "Mathematical 
Expectation; General Random Variables") 

The pair {X, Y} has joint density fxY {t, u) = 1/2 on the square region bounded hy u = 1 + t, 
u = 1 — t, u = 3 — t, and u = t— 1. 

W={ ] ~ =Iq{X,Y)X + Iqc{X,Y)2Y 

2Y formax{A:, Y} > 1 

where Q = {{t, u) : max{t, «}<!} = {{t, u) : t <l,u <1}. 
Determine E[W] and Var [W]. 
ANALYTIC SOLUTION 

The intersection of the region Q and the square is the set for which Q <t <1 and l — t<u<l. 
Reference to Figure 11.3.2 shows three regions of integration. 

E[W] = - / tdudt+- / 2ududt+- / 2u dudt = 11/6 ^ 1.S333 
2 Jo Ji-t 2 Jo Ji 2 Ji Jt-i 

E[W^]=- / t^dudt+- / 4u^dudt+- / 4u^ dudt = 103/24: 
2 Jo Ji-t ^ Jo Ji 2 Ji Jt-i 

Var [W] = 103/24 - (11/6)^ = 67/72 w 0.9306 



^^"Mathematical Expectation; General Random Variables", Example 15: A function with a compound definition 
<:http://cnx.org/content/m23412/latest/#fs-id7964278> 



http://cnx.Org/content/m23441/l.7/ 



Connexions module: m23441 



tuappr 

Enter matrix [a b] of X-range endpoints [0 2] 

Enter matrix [c d] of Y-range endpoints [0 2] 

Enter number of X approximation points 200 

Enter number of Y approximation points 200 

Enter expression for joint density ((u<=min(t+l,3-t))& 

(u$gt;=max(l-t,t-l)))/2 
Use array operations on X, Y, PX, PY, t, u, and P 
M = max(t,u)<=l; 
G = t.*M + 2*u.+(l - M) ; 
EG = total(G.*P) 
EG = 1.8340 

VG = total(G.~2.*P) - EG~2 
VG = 0.9368 
[Z,PZ] = csort(G,P) ; 
EZ = Z*PZ' 
EZ = 1.8340 
VZ = (Z.~2)*PZ' - EZ~2 
VZ = 0.9368 



Z = g(X,Y) 
E[g(X,Y)] 
'/. Theoretical 11/6 = 1.8333 

Theoretical 67/72 = 0.9306 
°/o Distribution for Z 
°/o E[Z] from distribution 



Example 8: A function with compound definition 

fxY (i, m) = 3 on < u < < 1 



Z = Iq {X, Y)X + Iqc {X, Y) for Q = {{t, u):u + t<l} 
The value tg where the Hne u = 1 — t and the curve u = meet satisfies = 1 — to- 

E[Z]=3 t dudt + 3 t I dudt + 3 / dudt = - {5to - 2) 
Jo Jo Jto Jo Jto Ji-t 4 

For E [Z"^] replace t by in the integrands to get E [Z"^] = (25io - 1) /20. 

Using to = (a/5 - 1) /2 « 0.6180, we get Var [Z] = (2125to - 1309) /80 w 0.0540. 
APPROXIMATION 

'/ Theoretical values 
to = (sqrt(5) - l)/2 
to = 0.6180 
EZ = (3/4)*(5*t0 -2) 
EZ = 0.8176 
EZ2 = (25*t0 - l)/20 
EZ2 = 0.7225 

VZ = (2125+tO - 1309)/80 

VZ = 0.0540 

tuappr 

Enter matrix [a b] of X-range endpoints [0 1] 
Enter matrix [c d] of Y-range endpoints [0 1] 
Enter number of X approximation points 200 
Enter number of Y approximation points 200 
Enter expression for joint density 3*(u <= t.~2) 
Use array operations on X, Y, t, u, and P 
G = (t+u <= l).*t + (t+u > 1); 
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Theoretical = 0.8176 
'/. Theoretical = 0.0540 



EG = total(G.*P) 
EG = 0.8169 

VG = total(G.~2.*P) - EG~2 
VG = 0.0540 
[Z,PZ] = csort(G,P) ; 
EZ = Z*PZ' 
EZ = 0.8169 
VZ = (Z.~2)*PZ' - EZ~2 
VZ = 0.0540 

Standard deviation and the Chebyshev inequality 

In Example 5^^ from "Functions of a Random Variable," we show that if X N (^jj,, <t^) then Z = ^^^zR 
N (0, 1). Also, E[X]=ii and Var [X] = a^. Thus 



PI \^LjA<t] =P{\X-fi\ <to) = 2$(t)-l 



(40) 

For the normal distribution, the standard deviation a seems to be a natural measure of the variation away 
from the mean. 

For a general distribution with mean /x and variance a^, we have the 
Chebyshev inequality 



\X-fi\ 



> a 



<\ or P{\X-fx\>aa) < \ 



(41) 



In this general case, the standard deviation appears as a measure of the variation from the mean value. This 
inequality is useful in many theoretical applications as well as some practical ones. However, since it must 
hold for any distribution which has a variance, the bound is not a particularly tight. It may be instructive 
to compare the bound on the probability given by the Chebyshev inequality with the actual probability for 
the normal distribution. 

t = 1:0.5:3; 
p = 2*(1 - gaussian(0 , 1 ,t) ) ; 
c = ones(l,length(t)) ./(t.~2) ; 
r = c./p; 

h = [' t Chebyshev Prob Ratio']; 

m = [t ; c ; p ; r] ' ; 
disp (h) 

t Chebyshev Prob 



disp (m) 

1.0000 
1.5000 
2.0000 
2.5000 
3.0000 

— □ 



1.0000 
. 4444 
0.2500 
0.1600 
0.1111 



0.3173 
0.1336 
0.0455 
0.0124 
0.0027 



Ratio 

3.1515 
3.3263 
5.4945 
12.8831 
41.1554 



DERIVATION OF THE CHEBYSHEV INEQUALITY 

Let A = {\X~ii\> aa} = {{X - > a^a^}. Then aV^J^ < {X - nf 

Upon taking expectations of both sides and using monotonicity, we have 



aV^P (A) <E {X-n) 



(42) 



"Functions of a Random Variable", Example 5: The normal distribution and standardized normal distribution 
^http://cnx.org/content/m23329/latest/#fs-idll69352255176> 



http://cnx.Org/content/m23441/l.7/ 
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from which the Chebyshev inequality follows immediately. 

— □ 

We consider three concepts which are useful in many situations. 
Definition. A random variable X is centered iff E [X] = 0. 

X = X — fj, is always centered. (43) 
Definition. A random variable X is standardized iff E [X] = and Var [X] = 1. 

X* = — = — is standardized (44) 

a a 

Definition. A pair {X, y} of random variables is uncorrelated iff 

E [XY] -E[X]E [Y] = (45) 

It is always possible to derive an uncorrelated pair as a function of a pair {X, Y}, both of which have finite 
variances. Consider 

U={X* + Y*) V={X* - Y*) , where X* = ^ ~ '^^ , Y* = ^-^J^ (46) 

ax cry 

Now E[U] = E [V] = and 

E [UV] =E{X* + Y*) {X* -Y*) = E Ux*)^] - E \{Y*f] =1-1 = (47) 

so the pair is uncorrelated. 

Example 9: Determining an uncorrelated pair 

Wo uso the distribution for Examples Example 10^^ from "Matliomatical Expectation: Simple 
Random Variables" and Example S {Z = g {X, Y) (Example 10^'' from "Mathematical Expectation: 
Simple Random Variables")), for which 

E [XY] -E[X]E [Y] ^ (48) 



jdemol 
jcalc 

Enter JOINT PROBABILITIES (as on the plane) P 
Enter row matrix of VALUES of X X 
Enter row matrix of VALUES of Y Y 

Use array operations on matrices X, Y, PX, PY, t, u, and P 
EX = total (t.*P) 
EX = 0.6420 
EY = total (u.+P) 
EY = 0.0783 
EXY = total (t . *u. *P) 
EXY = -0.1130 
c = EXY - EX+EY 

c = -0.1633 '/. {X,Y} not uncorrelated 



^^"Mathematical Expectation: Simple Random Variables", Example 10: Expectation for Z = g{X, Y) 
<http://cnx.Org/V.ontent/m23387/latest/#fs-idll64693552188> 

"Mathematical Expectation: Simple Random Variables", Example 10: Expectation for Z = g{X, Y) 
<http://cnx.Org/content/m23387/latest/#fs-idll64693552188> 
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VX = total (t.~2.*P) - EX~2 
VX = 3.3016 

VY = total(u.~2.*P) - EY~2 

VY = 3.6566 

SX = sqrt(VX) 

SX = 1.8170 

SY = sqrt(VY) 

SY = 1.9122 

X = (t - EX)/SX; °/o Standardized random variables 

y = (u - EY)/SY; 

uu = X + y; y, Uncorrelated random variables 

vv = X - y; 

EUV = total (uu.+vv.*P) % Check for uncorrelated condition 

EUV = 9.9755e-06 Differs from zero because of roundoff 



http://cnx.Org/content/m23441/l.7/ 



